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Abstract 

The major focus of this article is to analyze the forced convective heat transfer in a steady boundary layer flow of Sisko fluid 
over a nonlinear stretching sheet. Two cases are studied, namely (i) the sheet with variable temperature (PST case) and (ii) 
the sheet with variable heat flux (PHF case). The heat transfer aspects are investigated for both integer and non-integer 
values of the power-law index. The governing partial differential equations are reduced to a system of nonlinear ordinary 
differential equations using appropriate similarity variables and solved numerically. The numerical results are obtained by 
the shooting method using adaptive Runge Kutta method with Broyden's method in the domain[0,oo). The numerical 
results for the temperature field are found to be strongly dependent upon the power-law index, stretching parameter, wall 
temperature parameter, material parameter of the Sisko fluid and Prandtl number. In addition, the local Nusselt number 
versus wall temperature parameter is also graphed and tabulated for different values of pertaining parameters. Further, 
numerical results are validated by comparison with exact solutions as well as previously published results in the literature. 
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Introduction 

The flow and heat transfer over stretching surfaces are relevant 
to many technological processes. The study of boundary layer flow 
and heat transfer induced by stretching heated surface has 
attracted considerable attention of many researchers. Some 
exemplary applications of such study are glass fibre production, 
the cooling and drying of paper while paper production, drawing 
of plastic films and so forth. Starting from the pioneering work of 
Sakiadis [1] numerous aspects of momentum and heat transfer 
over a stretching sheet have been considered. Some investigations 
dealing with the flow and heat transfer over nonlinear stretching 
surface are reported in Refs. [2-6]. 

All of the above mentioned studies are restricted to viscous 
fluids. However, there is a voluminous body of knowledge that 
testifies the wide occurrence of non-Newtonian fluids in industrial 
sector including food, pharmaceutical, polymer and plastic, 
mineral suspensions, cosmetics, personal care products, toiletries, 
construction materials and biological products [7]. In view of 
overwhelming practical utility of non-Newtonian fluids, several 
researchers have proposed numerous models for non-Newtonian 
fluids, including the Sisko fluid model [8]. The Sisko fluid is of 
much significance due to its adequate description of many non- 
Newtonian fluids over the most important range of shear rates. 
The three parametric Sisko model can be considered as a 
generalization of Newtonian and power-law fluids. The three 
constants in the model can be chosen with great ease for specific 
fluids and model is found to be good in predicting the shear 
thinning and shear thickening behaviors. 

Extensive research has been commenced over past few years on 
the flow and heat transfer of non-Newtonian fluids over a 



stretching surface owing to its tremendous industrial utilization. 
Firstly, Schowalter [9] obtained the similar solutions for the 
boundary layer flow for power-law pseudoplastic fluids. Jadhav 
and Waghmode [10] analyzed the heat transfer to power-law fluid 
over a permeable flat plate with heat flux boundary condition. 
Howel et al. [11] considered the laminar flow and heat transfer of a 
power-law fluid over a stretching sheet. Hassanien et al. [12] 
numerically analyzed the flow and heat transfer for a stretching 
sheet for non-uniform temperature distribution. Able et al. [13] 
studied the flow and heat transfer to a power-law fluid over a 
stretching sheet by considering variable thermal conductivity and 
heat source. 

Polymeric suspensions such as waterborne coatings are identi- 
fied to be non-Newtonian in nature and are proven to follow the 
Sisko fluid model [14]. The Sisko fluid model was originally 
proposed for high shear rate measurements on lubricating greases 
[15]. Khan et al. [16] examined the steady flow and heat transfer 
of a Sisko fluid in annular pipe. Then, Khan and Shahzad [17,18] 
developed the boundary layer equations for Sisko fluid over planer 
and radially stretching sheets and found the analytical solutions for 
only integral values of the power-law index. The utmost studies 
relating to the heat transfer of Sisko fluid involve only one 
dimensional flows and literature survey indicates that no work has 
so far been communicated with regards to heat transfer in a 
boundary layer flow for Sisko fluid over a nonlinear stretching 
sheet with variable surface temperature and variable heat flux. 

The prescribed surface temperature and prescribed surface heat 
flux are the generalizations of constant surface and constant heat 
flux thermal boundary conditions. Such boundary conditions are 
encountered in polymer processing industry, where the surface 
temperature may be an arbitrary function of time or space. The 
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prescribed heat flux may be taken when there is surface heat 
generation via solid-solid friction, as in frictional welding and 
conveying of solids in screw extruders, is an example. Moreover, 
certain types of intensive radiation or convective heating that are 
weak functions of surface temperature can also be treated as a 
prescribed surface heat-flux boundary condition [19]. 

This paper brings to focus the heat transfer in a boundary layer 
flow of Sisko fluid over a nonlinear stretching sheet with variable 
wall temperature and heat flux boundary conditions. Khan and 
Shahzad [18] in their work only considered the integral values of 
the power-law index. In this investigation, the non-integer values 
are also taken into consideration. 

Formulation of the Problem 

Flow Analysis 

We consider the steady, laminar and incompressible flow of a 
Sisko fluid impinging normal to a stretching sheet coinciding with 
the plane y = 0. The flow is induced due to stretching of the sheet 
along the x-direction while keeping the origin fixed with velocity 
U w = cx s , where cand s are positive real numbers relating to 
stretching of sheet. For a two-dimensional flow, we assume the 
velocity field of the form 

V=[u(x,y),v(x,y),0], (1) 

Where (x,y) denotes the Cartesian coordinates along and 
perpendicular to the sheet, u and v are the velocity components of 
fluid along the x — and y — directions, respectively. 

Using the simplifying assumptions, the momentum equation 
characterizing the steady boundary layer flow takes the form 
(see ref. [18] for details) 

/ du du\ d 2 u d ( du\ n 
Where a, b and «(>0) are the material constants. 



The flow is subject to the boundary conditions 

U = cx s ,v = 0 at y = 0, (3) 



w->0 as >>->oo. (4) 
In view of the similarity transformations [18] 

rj=^Re^ and ^(x,y)=xURe^~ l f(l), (5) 

and after mathematical simplification, we obtain the following 
problem 

Ar+n[-ff- 1 r+ ( * (2 7 + \ )+1 )//"-*(n 2 =o, (6) 

/(0)=0,/'(0) = l,/'(a>) = 0, (7) 

where prime denotes the differentiation with respect to f/and the 
dimensionless quantities are defined by 

2 

Re a = pxU/a, Re b = px n U 2 ~ n /b and A = Re n b + l / Re a . (8) 

The physical quantity of major interest is the local skin friction 
coefficient and in dimensionless form is given by [18] 

l -Re^C f = Af"(0)-[-f"(0)] n . (9) 



Heat Transfer Analysis 

By using the usual thermal boundary layer approximations, 
neglecting viscous dissipation and heat generation, the energy 
equation for temperature field T— T(x,y) is given by 




Figure 1. The temperature profile 9{r\) in the PST case for different values of the power-law index n when s=\.5, Pr =3,m=l 
and A = 1 are fixed. 

doi:1 0.1 371 /journal.pone.01 00056.g001 
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Figure 2. The temperature profile g(rj) in the PHF case for different values of the power-law index n when 5=1.5, Pr =3, m=l 
and A = 1 are fixed. 

doi:1 0.1 371 /journal.pone.01 00056.g002 

k 

where a = is the thermal diffusivity with c p as the specific heat 

pc P 

of fluid at constant pressure and K the thermal conductivity. 

Two kinds of thermal boundary conditions are considered in 
heat transfer analysis, and they are treated separately in the 
following sections. 



dT dT d 2 T 

U h V = (X 

dx dy dy 2 



(10) 
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Figure 4. The temperature profile g(rj) in the PHF case for different values of the stretching parameter 5 when m= 1.5, Pr =3 and 
A = 1 are fixed. 

doi:1 0.1 371 /journal.pone.01 00056.g004 




Figure 5. The temperature profile 9(r]) in the PST case for different values of the material parameter A when 5= 1.5 and m= 1 are 
fixed. 

doi:1 0.1 371 /journal.pone.01 00056.g005 



PLOS ONE | www.plosone.org 



4 



June 2014 | Volume 9 | Issue 6 | e1 00056 



Convective Heat Transfer in Sisko Fluid 



2* 



(■) 

n-0.5 




A- 0.0 

A- 1.0 

A- 2.0 

A - 3,0 



■ i i I 




A -0.0 

A-1.0 

- A -2.0 
A- 3.0 




A -0.0 

A- 1.0 

A-2J0 

A -3.0 



Figure 6. The temperature profile g{t]) in the PHF case for different values of the material parameter A when s = 1.5 and m= 1 are 
fixed. 

doi:1 0.1 371 /journal.pone.01 00056.g006 




Figure 7. The temperature profile Q(rj) in the PST case for different values of the temperature parameter m when s= 1.5, Pr = 1 and 
A = 1 are fixed. 

doi:1 0.1 371 /journal.pone.01 00056.g007 
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Figure 8. The temperature profile g(r\) in the PHF case for different values of the temperature parameter m when s= 1.5, Pr = 1 and 
A = 1 are fixed. 

doi:1 0.1 371 /journal.pone.01 00056.g008 



Prescribed Surface Temperature (PST Case) 

In this case, the thermal boundary conditions are 



The boundary conditions for 6(rj) follow from Eqs. (11) and 
(12) are 



T=T ™(= T ™+ B (l) m ) at J = 0, (11) 



T^Too as j>->oo, (12) 

where B is a constant, L the characteristic length, the constant 
fluid temperature far away from the sheet, T w the wall 
temperature and m the wall temperature parameter. 
Defining the non-dimensional temperature 0(rj) by 



T-T a 



and using Eqs. (5), (11) and (12), Eq. (10) takes the form 



9"+Pr ^ (2w - 1 , )+1 /e'-^)=0, 



n + \ 



(13) 



(14) 



xUR 



2 



where Pr = — is the generalized Prandtl number. 



0(0) = 1, and %)^0 as ^oo. 

The rate of heat transfer at the sheet surface is 

q w \ y=0 = -K(T w -T x )9'(0)^Rl /n + l , 
which on simplification reduces to 

Re- ll " +x Nu x =-6'{0). 



where Nu x = 



xq w 



k(T w Tix) ) 



is the local Nusselt number. 



(15) 



(16) 



(17) 



Prescribed Heat Flux (PHF Case) 

In PHF case, the dimensionless temperature g(t]) is defined as 



g{n)-- 



T-T a 



(D/K)(x/LrxR^" +l ' 



(18) 



with the following boundary conditions 
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Figure 9. The temperature profile 6(rj) in the PST case for different values of the Prandtl number Pr when s=\.5,m= \ and A = \ are 
fixed. 

doi:1 0.1 371 /journal.pone.01 00056.g009 



q w =-K 



dT 



y = 0 



T^-T^asy^ oo, 



(19) 



(20) 



where D is a constant, and m = 0 leads to a constant heat flux case. 
Using Eqs.(5) and (18) in Eq. (10), one finds that 



gt!+ Pr f^^^ + fc^ /V(M+1)/ g] = 0 , (21) 



ft+1 



ft+1 



g'(0) = — 1 and g(rj)^>Q as rj^>co. 
The local heat transfer coefficient is 



h~- 



(22) 



(23) 



resulting in the local Nusselt number given by 

1 



Rl +X Nu x = 



(24) 



hx 



with Nu x = — as the local Nusselt Number. 
k 



Exact Solutions for Particular Cases 

Since Eq. (6) has simple exact solution to a special case, namely 
n = s= 1 (see ref. [18] for detail). For this case with m = 0, Eq.(14) 
reduces to 



(25) 



The exact solution of Eq. (25), in terms of the incomplete 
Gamma function, satisfying boundary conditions (15) is 



%) = 



Pr Pr Pr 

y y y 



,.Pr 



Pr Pr 



(26) 



r(-,,o)-r(-,,-,) 

y y y 
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Figure 1 1. A comparison of the exact and numerical results (solid line numerical results and open circles exact results) in the PST 
case when n= \, s = 1 and A = 1 are fixed. 

doi:1 0.1 371 /journal.pone.01 00056.g01 1 



PLOS ONE | www.plosone.org 



8 



June 2014 | Volume 9 | Issue 6 | e1 00056 



Convective Heat Transfer in Sisko Fluid 




Figure 12. A comparison of the exact and numerical results (solid line numerical results and open circles exact results) in the PST 
when n= \, s = land A=() are fixed. 

doi:1 0.1 371 /journal.pone.01 00056.g01 2 

where /?= = and r( ) the incomplete Gamma function. 

x _ ft M[Pr-/w, Pr + ^-Pre-*] 
For n = s= 1 and A = 0,Eq. (14) reduces to Svl) — Yx e M[Pr —m Pr + 1 — Pr] ' 



O"+Pr(fO , -mf0) = O. 



(27) 



The exact solution of the above equation satisfying the 
boundary conditions (15), in terms of Kummer's function, is 
expressed as [20] 



W M[Pr-m,Pr+l,-Pr] ' 1 } 



Where M( ) is the Kummer's function. 

When n = s= 1 and A = 0, Eq.(21) has exact solution satisfying 
boundary conditions (22) of form [20] 



Solution Methodology 

In general it is very difficult to find the exact analytical solution 
of non-linear two point boundary value problem (6), (14) and (21) 
along with boundary conditions (7), (15) and (22). Therefore, these 
problems are solved numerically by the shooting technique. The 
equations are firstly written as a system of five first order ordinary 
differential equations. Then the corresponding initial value 
problems are solved by the adaptive Runge-Kutta method. The 
initially guessed values f"(0) and 0'(O) are refined iteratively 
using the Boryden's method to satisfy boundary condition at 
infinity. The iterative process is terminated when the absolute 
error is less than the tolerance 10 -4 . 




Figure 13. A comparison of the exact and numerical results (solid line numerical results and open circles exact results) in PHF case 
when n = 1, s = 1 and A = 0 are fixed. 

doi:1 0.1 371 /journal.pone.01 00056.g01 3 
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Figure 14. The variation of the local Nusselt number with m for some value of the power-law index n and Prandtl number Pr. 

doi:10.1371/journal.pone.0100056.g014 



Numerical Results and Discussion 

The main aim of the present study is to investigate the heat 
transfer to a Sisko fluid over a nonlinear stretching sheet with the 
non-isothermal wall temperature and variable heat flux boundary 
conditions. In order to comprehend the heat transfer phenomena 
in detail, Eqs. (14) and (21), with their respective boundary 
conditions, are solved numerically and the results are displayed 
through graphs. Mainly, the effects of the power-law index n, 
stretching parameters, Prandtl number Pr, material parameter^, 
and surface temperature parameter m are investigated in detail 
both for PST and PHF cases. Moreover, the heat transfer aspects 
are explored in terms of the local Nusselt number at the wall. 
Figures 1-10 give the various perspectives of heat transfer for both 
PST and PHF cases. The values of the local Nusselt number are 
recorded in tables 1 and 2 for both PST and PHF cases, 
respectively. 

The effect of the power-law index n on temperature profile 6(rj) 
is presented in figures 1 and 2 for both PST and PHF cases, 
respectively. A wider range of the power law index n is spanned in 
the analysis (0.3 <n<2.5) covering from highly shear thinning to 
thickening fluids. 

Most polymeric melts and solutions have value of the power law 
index in the range 0.3-0.7 [21], moreover, typical starch-in-water 
and corn flour-in-water dispersion have values of power-law index, 
1.4-1.6 [22]. The figure \{a) shows that the temperature profile 
increases when the power-law index is decremented below one. 
This increase can be attributed to the fact that the reduction in 



viscosity favors the heat transfer. The increase is large for smaller 
values of n. The temperature profile also decreases as the value of 
n is incremented gradually (figure \b) for n> 1. A comparison of 
the two patterns reveals a strong dependence of heat transfer on 
flow behavior index for n < 1 and weaker when n > 1 . The figures 
2(a, b) elucidate the temperature profiles for the PHF case for 
n < 1 and n > 1 , respectively. It is quite clear from these figures that 
the temperature profiles decrease when the value of the power-law 
index is increased for shear thinning (n < 1 ) and shear thickening 
(n>\) regimes. These figures show that, in the PHF case, the 
behavior of temperature profile with change in n is similar to the 
PST case. However, a stronger dependence of the temperature 
profile can be noticed in the PHF case. Both the figures put in an 
evidence that an augmentation in the value of the power-law 
results in a decrease of the thermal boundary layer thickness. 

The stretching parameter s affects the temperature profile and 
thermal boundary layer by virtue of imparting shear stress at the 
boundary. A broader range 0.5<s<2.0 is included in the 
analysis, including s = 1.0 for linear stretching. Its effects on heat 
transfer for both boundary conditions are presented in figures 3 
and 4 respectively. For the PST case, when the value of the power- 
law index (n< 1), the effect of the stretching parameter s on heat 
transfer to the fluid is very meager; however, the thermal 
boundary layer thickness increases a bit as the value of s is 
incremented progressively. We can ascribe this increase to the 
amplification of shear stress with each increment in s. For shear 
thinning regime (n< 1) the enhancement of shear stress lowers the 
effective viscosity, that favors the heat transfer. Although, for n = 1 , 



Table 1. The numerical values of the Local Nusselt number in the PST case when m = 1 and Pr = 1 are fixed. 





-Re b 


l/n+l Nu x 








A 


5 


A7=0.5 


A7= 1.0 


A7= 1.5 


1 


1 


1.01162 


1 .07693 


1.10805 


2 


1 


1.04162 


1.11064 


1.14661 


3 


1 


1.05911 


1.12991 


1.16897 




2 


0.96658 


1.14410 


1.23540 


2 


2 


1 .00436 


1.18608 


1.28120 


3 


2 


1.02653 


1.21105 


1.30854 



doi:1 0.1 371/journal.pone.01 00056.t001 
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Table 2. The numerical values of the Local Nusselt number in the PHF case when m = 


1 and Pr = 1 are fixed. 




Re; l/n+l Nu x 


A 5 


A7=0.5 


77=1.0 


77=1.5 


1 1 


0.86916 


1.07706 


1.17826 


2 1 


0.89766 


1.11033 


1.21713 


3 1 


0.91446 


1.13041 


1.23968 


1 2 


0.27308 


0.97012 


1.23540 


2 2 


0.29081 


1.00985 


1.28121 


3 2 


0.30209 


1.03362 


1.30854 





doi:1 0.1 371/journal.pone.01 00056.t002 



the effect of s on heat transfer is noticeable (figure 3/3), but a 
contrary behavior for the temperature and thermal boundary layer 
thickness is observed. Figure 3(c) presents the effect of the 
stretching parameter s on the temperature profile for n= 1.5. The 
temperature and thermal boundary layer decrease and s affects the 
thermal boundary layer thickness substantially for this case. Figure 
4(a) depicts that the temperature profile raises as the value of s is 
augmented for shear thinning (n<\) regime. The temperature 
near the wall grows quite prominently as the value of s is 
incremented. The character of temperature profile, when 
n=l, and stretching parameter s varied is shown in figure 4(b). 
This figure reveals that on increasing the value of the stretching 
parameter s, the temperature profile increases near the wall and 
then decreases away from the wall. For shear thickening (n= 1.5) 
regime, the temperature profile and corresponding thermal 
boundary layer show a diminution as the value of s is raised 
(figure 4c). 

The effect of the material parameter A on the temperature 
profile for nonlinear stretching is presented in figures 5 and 6 
respectively, for PST and PHF cases. These figures also make a 
comparison amongst the temperature profiles of the Newtonian 
fluid 04 = 0 and n = 1) and the power-law fluid (^4 = 0 and n ^ 1) 
with those of the Sisko fluid (^4^0). Figures 5(a — c) reveal that 
the temperature profile and corresponding thermal boundary 
layer thickness depress in each case with increasing value of the 
material parameter A. It can be noticed from the sketches that the 
effect of A on the thermal boundary layer is pronounced for low 
Prandtl number as compare to higher one. Figures 6(a — c) 
demonstrate that the temperature profile and corresponding 
thermal boundary layer thickness also decrease with an increase 
in the value of A for different values of power-law index n. Effects 
are more prominent for lower Prandtl numbers. 

The influence of the wall temperature parameter m for both 
PST and PHF cases on the temperature distribution and thermal 
boundary layer thickness is presented in figures 7 and 8 



respectively. Numerical solutions are sought in the range 
-0.5<ra< 1.0 and 0<ra< 1.5 for PST and PHF cases, respec- 
tively. These ranges are selected on the basis of physically 
acceptable solutions. From these figures, it is observed that as the 
value of m is incremented progressively from negative to positive, 
the temperature and thermal boundary layer thickness decrease as 
shown for some values of the power-law index n But, this effect 
turn out to be diminish for larger values of n. 

The Prandtl number Pr of a fluid plays a dominant role in 
forced convective heat transfer. The computations are carried out 
for large Prandtl number, many non-Newtonian fluids exhibit a 
value of Prandtl number as high as 100 or even greater [23,24]. Its 
effect for the PST case for different values of power-law index n is 
shown in figures 9(a,b,c). It is notified that the heat transfer 
process is augmented prominently by thinning the thermal 
boundary layer thickness whenPris increased. The augmentation 
can be ascribed to the enhanced momentum diffusivity for larger 
Prandtl number and the heat transfer mainly occurs due to 
advection. The same qualitative aspects are observed for the PHF 
case as shown in figures \0(a,b,c). 

The numerical results are also compared with those of exact 
ones as special cases of the problem (figures 11-13). An excellent 
agreement confirms the credibility of our numerical solutions. 

Figure 1 4(a) depicts the variation of the local Nusselt number 
6'(0) with change in wall temperature parameter m. It is clearly 
observed from the plots that the heat transfer at the wall increases 
rapidly for the fluids with greater Pr. Moreover, the fluids with 
n > 1 carry larger heat out of heated surface. These plots also give 
the value of mwhere the value of 6'(0) approaches to zero. This 
value of m is strong function of the power-law index n. Figure 
14(7?) illustrates the effect of m on the local Nusselt number for the 
Sisko fluid (both for shear thinning and shear thickening regimes). 
The Prandtl number affects the heat transfer in same fashion as 
previous one, but here more heat is transferred for each value of m 
as compared to the PST case. 



Table 3. A comparison of the local Nusselt number in the PHF case. 





n 


Chen [25] 




Present study 






Re- x,n+l Nu x 


g(0) 


Re- l/n+l Nu x 


0.5 


0.659092 


1.517239 


0.659091 


1.517019 


1.0 


0.590316 


1 .694009 


0.590312 


1.694001 


1.5 


0.561785 


1 .780041 


0.561781 


1.780040 



doi:1 0.1 371/journal.pone.01 00056.t003 
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Table 1 briefly shows the trend in the local Nusselt number for 
shear thinning and thickening fluids when the material parameter 
A and stretching parameter s are varied for the PST case. It can be 
noticed from this table that the local Nusselt number shows a boost 
with an each increment of A. It is also clear that there is much 
boost in the local Nusselt number for flow behavior index n > 1 . 
Table 2 displays the same informations and qualitatively similar 
trends are observed for the PHF case. Also, in table 3, a 
comparison with the obtained results of Chen [25] is made and 
excellent agreement is found. 

Conclusions 

The analysis in this article was focused on heat transfer in a 
Sisko fluid over a nonlinear stretching sheet both for variable 
temperature and heat flux boundary conditions. Numerical 
solutions are obtained and verified as special cases of the problem 
with exact results. Our computations have indicated that: 

• A rapid decrease in the thermal boundary layer thickness was 
observed for Sisko fluid with shear thinning properties with 
PHF boundary condition. 

• Thermal boundary layer thickness was increased for power- 
law-index n < 1 , and decreased for n > 1 when stretching 
parameter swas increased. The effects were more noticeable 
for PHF case. 
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